Time evolution of averages 
in dynamical systems driven by noise 



> 



Andrey L. Pankratov 

Institute for Physics of Microstructures of RAS, GSP 105, 
Nizhny Novgorod, 603600, RUSSIA. E-mail: alp@ipm.sci-nnov.ru 

This paper presents description of time evolution of averages of Markov process in wide range 
of noise intensity. Exact expression of time scale of average evolution has been obtained. It has 
been demonstrated numerically that for purely noise-induced transitions (transitions over potential 
barriers) the time evolution of mean coordinate is a simple exponent with a good precision even for 
the case when the potential barrier height is comparable or smaller than the noise intensity. Also it 
has been demonstrated that nonlinear system may be "linearized" by a strong noise. 
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Overdamped Brownian motion in a field of force (Markov process with nonlinear drift coefficient) 
is a model widely used for description of noise-induced transitions in different polystable systems. 



I. INTRODUCTION 



Values that may be observed in experiment are different averages and knowledge of their time 
evolution is a subject of great theoretical and practical importance in a wide variety of tasks (e.g. 
tasks of Josephson electronics |J, stochastic resonance 0, ratchet effect @|j 0> an d so on). 

In paper by Nadler and Schulten || the authors have demonstrated exponential behavior of ob- 
servables for relaxation processes in systems having steady states. Namely, using the approach of 
"generalized moment approximation" the authors of || obtained some exact characteristic time 
scales and for particular example of a rectangular potential well demonstrated good coincidence of 
kQ| exponential approximation with numerically obtained observables. 

In the frame of this paper we consider time evolution of averages of an overdamped Brownian 
motion. Extending the approach by Malakhov we have obtained exact characteristic scales of time 
evolution of averages (valid for arbitrary potentials and noise intensity). Some particular examples 
of time evolution of mean coordinate of overdamped Brownian motion in monostable and bistable 
potentials have been considered. It has been demonstrated that time behavior of mean coordinate 
is very well approximated by the exponent for purely noise-induced transitions (transitions over 
potential barriers) even in the case of a large noise intensity in comparison with potential barrier 
height. Also, an effect of "linearization" of nonlinear systems by a strong noise has been observed. 

It is necessary to mention that in difference from approach used in ||, our approach allows to 
describe time evolution of averages for tasks with arbitrary boundary conditions, e.g. the case of 
periodic boundary conditions (p(n) = (p(—n) as application to tasks of Josephson electronics will 
be soon presented elsewhere. Moreover, as it seems now, the use of the presented approach for 
description of ratchet effect may allow to get the required optimal frequency of flashing H for 
periodic potentials of arbitrary shape. 

II. MAIN EQUATIONS AND SET UP OF THE PROBLEM 

Consider a process of Brownian diffusion in a potential profile <P(x). Let a coordinate x(t) of the 
Brownian particle described by the probability density W(x,t) at initial instant of time has a fixed 
value x(0) = x , i.e. the initial probability density is the delta function: 
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W(x, 0) = 6(x-x ). 



(1) 



It is known that the probability density W(x, t) of the Brownian particle in the overdamped limit 
(Markov process) satisfies to the Fokker-Planck equation (FPE): 

dW{x,t) dG(x,t) _ 1 f d 
dt dx B 1 dx 

Here B = h/kT, G(x, t) is the probability current, h is the viscosity (in computer simulations we put 
h — 1), T is the temperature, k is the Boltzmann constant and <p(x) = <P(x)/kT is the dimensionless 
potential profile. In this paper we will only consider an overdamped Brownian motion in potentials 
ip(x) such that </?(±oo) = +oo, i.e. eventually the probability density will reach a steady-state 

W s t(x) = Ae~ v ^ x ' ^ (A = 1/ / e~ v ^dx). This leads to the following boundary conditions: 

— oo 

Gr(±oo, t) = 0. Besides, the growth of the walls of the potential <p(x) should be fast enough for 
applicability of the method described below (for (p(x) ~ x n , where n < 1, the method not always 
gives correct results). 
We will search for the average rrif(t) in the form: 

+oo 

m f (t) =< f(x) >= J f(x)W(x,t)dx. (3) 

— oo 

Let us mention that the approach we use may be easily extended for description of time evolution 
of some kind of conditional averages when the limits ±oo may be substituted by arbitrary interval 
(c, d), such that boundary conditions at c and d are arbitrary. 

Following this obvious definition (^) it is necessary to know the solution of the FPE - the probability 
density W(x,t) - for obtaining the required average. But the nonstationary solution of FPE is 
unknown in general case of nonlinear systems. However, as it has been demonstrated recently in 

d 

H, |7j], the time evolution of the probability P(t) = J W(x,t)dx may be very well (with only few 

c 

percent mistake) approximated by the exponent even when the noise intensity is larger than a 
potential barrier height. Besides, it is well-known that the time evolution of mean coordinate of 
Brownian motion in a parabolic potential (linear system) is exactly exponential. This allows to hope 
that time evolution of different averages of Markov process may be in many cases very well described 
by exponential approximation in wide range of noise intensity. 

Modifying the approach by Malakhov ||, 0, it is possible to define some characteristic scale of time 
evolution of the required average (similar way was used in ||) and obtain exact analytic expression 
for this time scale. After that, substituting the obtained time scale into the factor of exponent and 
comparing the obtained time evolution of the average with computer simulation results we will test 
applicability of this approximation for different types of potentials in wide range of noise intensity. 



dtp(x) 
dx 



W(x,t) 



+ 



d 2 W(x,t) 
dx 2 



(2) 



III. CHARACTERISTIC TIME SCALE OF EVOLUTION OF AVERAGE OF MARKOV PROCESS 

Let us consider attentively the FPE @. This is a continuity equation: 

dW(x,t) _dG(x,t) 
dt dx 



(4) 



To obtain necessary average rrif(t) @ let us multiply equation (|J) by the function f(x) and integrate 
it with respect to x from — oo to +oo. Then we get: 
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drrif(t) 
dt 



f(x)dG(x,t). (5) 



This is already ordinary differential equation of the first order (f(x) is known deterministic function), 
but nobody knows how to find G(x,t). 

Let us define the characteristic scale of time evolution of the average rrif(t) by similar way as it 
was done for time scale of the probability evolution in ||, |J: 

oo 

/ [mf(t) — m/(oo)] dt 
T ' M = ° m ,(0)- m/ (oo) ' (6) 

Definition (|]) is general in the sense that it is valid for any initial condition. But in the frame of 
this paper we restrict ourselves by the delta-shaped initial distribution ([I]) and consider 17 (x ) as 
function of xq. For arbitrary initial distribution the required time scale may be obtained from t/(xo) 
by simple averaging over initial distribution. If the required function m/(t) evolves exponentially in 
time Trifit) ~ e~'/ r then the time scale r in the factor of exponent coincides with tj{xq) defined by 

(!)• 

It is necessary to mention that definition (^j) gives correct results only for monotonically evolving 
functions rrif(t). Besides, rrif(t) should fast enough approach its steady-state value m^(oo) for 
convergence of the integral in (^). 

The required evolution time of the average rrif(t) may be obtained via the approach proposed 

by Malakhov ||, ||. This approach is based on the Laplace transformation method rrif(s) = 

00 

/ rrif(t)e~ st dt. We will only slightly modify this approach with respect to our task. 


Performing Laplace transform of formula (^), Eq. (Laplace transformation of gives: 

+00 

srrif(s) — mf(p) = — / f(x)dG(x, s)) and combining them, we get: 

—00 

+00 

z \ , x mt (0) — nit (00) — f f(x)dG(x,s) 

/ x ,. smAs) -m/(oo) ,. /w /v 7 -e» 

r*(x ) = lim . *) { ^ — (r = Inn = ^- — . (7) 

A u; s^o s [m f (0) -m/(oo)] »-o s[m/(0) -m/(oo)] w 

where s) is the Laplace transformation of the probability current G(x, s) — J G(x,t)e~ st dt. 



Following the approach by Malakhov, one can introduce the function H(x,s) = sG(x,s), and 
expand it in the power series in s: 

H(x, s) = sG(x, s) = H Q (x) + sHi(x) + s 2 H 2 (x) + ... (8) 

It is possible to find the differential equations for H n (x) (see 0): 

dHl(x) txt / \ e/ \ 

, = -W st (x) + 5(x - x ), 

dx (g) 
d 2 H n (x) dipjx) dH n (x) M n -o*A 

— ~T~^ 1 3 j - iJiln-i I , n = 2,6,4, ... 

dx z dx dx 

+00 

Here W st = Ae~^ x \ A = 1/ J e _I/P ^^(ix. Using the natural boundary conditions G(±oo,t) = 0, one 

—00 

X 

can obtain from (||]) H±(x) = —A J e~ v ^"'dv + l(x — Xq) and 
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X 



Ho(x) = B 



e-^ v) J e v>iy) [l{y-x )-F{y)]dydv- 

— oo 

-F(x) J e -<"« J e*M[l{y - x ) - F(y)]dydul , (10) 

— oo — oo ) 

where 

M +00 

F{u) = f e-^ v) dvj f e'^dv. (11) 



Here we restricted ourselves by obtaining of H 2 (x) only, because substituting the set @ into 
formula ([!]) one can get: 

+00 

/ f(x)dH 2 (x) 

TfM = ~ K(0)-m /( oo)] - (12) 

Thus, substituting H 2 (x) fllCf ) into formula (^2j) one can get the characteristic scale of time evolu- 
tion of any average rrif(t) for arbitrary potential such that <^(±oo) = 00: 

TfM = m f (0) -m f (oo) X 



x { I f(x)e-^ x) J e^ u) F(u)dudx -A J f{x)e~^ x) dx J e~^ x) J e^ u) F(u)dudx+ (13) 

X3 XQ —OO —OO 20 

OO OO X OO X 

+A f f{x)e- v(x) dx f e~^ x) f e^dudx - f f(x)e~^ x) f e^ u) dudx 



— OO XQ Xq Xq Xq 

+00 

where F(x) is expressed by ( pi]) and A = 1/ J e~^ x ^dx. 

—00 

Once we know the required time scale of the evolution of average, we can present the required 
average in the form: 

m f {t) = (to/(0) - m/(oo)) exp(-t/r/(a;o)) + m/(oo). (14) 

The applicability of this formula for several examples of time evolution of mean coordinate m(t) =< 
x(t) > will be checked in the next section. 

IV. TIME EVOLUTION OF MEAN COORDINATE OF MARKOV PROCESS 

As an example of description presented above, let us consider time evolution of mean coordinate 
of Markov process: 

+00 

m(t)—<x(t)>— / xW(x,t)dx. (15) 
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The characteristic time scale of evolution of the mean coordinate in general case may be easily ob- 
tained from (|T3"D substituting x for f(x). But for symmetric potentials <p(x) = <^{—x) the expression 
of time scale of mean coordinate evolution may be significantly simplified (to(oo) = 0): 



r m (x ) = — 




f^dudx 



(16) 



If x = then, as not difficult to check, T m (x ) = 0. 

Let us consider now some examples of symmetric potentials and check applicability of exponential 
approximation: 



First we should consider time evolution of mean coordinate in monostable parabolic potential 
f(x) = ax 2 /2 (linear system), because for this case time evolution of mean is known: 



where a = a'/kT and B = h/kT, so r par = B/a for linear system and does not depend on noise 
intensity and the coordinate of initial distribution xq. On the other hand, r m (xo) is expressed 
by formula fll6|) . Substituting parabolic potential <f(x) = ax 2 /2 in formula (|16"D, making simple 
substitutions and changing order of integrals, it can be easily demonstrated that r m (x ) = B/a = 
h/a' = T par , so for purely parabolic potential the time of mean evolution (|T6|) is independent of both 
noise intensity and xo as it must. This fact proves the correctness of the used approach. 

The second considered example is described by the monostable potential of the fourth order: 
(f(x) = ax 4 /4. In this nonlinear case the applicability of exponential approximation significantly 
depends on the location of initial distribution and the noise intensity. Nevertheless, the exponential 
approximation of time evolution of the mean gives qualitatively correct results and may be used as 
first estimation in wide range of noise intensity (see Fig. 1, a — 1). Moreover, if we will increase noise 
intensity further, we will see, that error of our approximation decreases and for kT = 50 we get that 
the exponential approximation and the results of computer simulation coincide (see Fig. 2, plotted 
in the logarithmic scale, a — 1, xq — 3). From this plot we can conclude that nonlinear system is 
"linearized" by strong noise - effect which is qualitatively obvious, but should be investigated further 
by the analysis of variance and higher cumulants. 

The third considered example is described by the bistable potential - the so-called "quartic" 
potential: <^(x) = ax A /4 — bx 2 /2. In this case the applicability of exponential approximation also 
significantly depends on the coordinate of initial distribution. If xq is far from the potential minimum, 
then there exists two characteristic time scales: fast dynamic transition to potential minimum and 
slow noise-induced escape over potential barrier. In this case exponential approximation gives not 
so adequate description of temporal dynamics of the mean, however may be used as first estimation. 
But if Xo coincide with the potential minimum, then the exponential approximation of the mean 
coordinate only few percent differs from results of computer simulation even in the case when noise 
intensity is significantly larger than the potential barrier height (strongly nonequilibrium case) (see 
Fig. 3, a = 1, b = 2, xq = 1.414). If, however, we consider the case when the initial distribution xq 
is far from the potential minimum, but the noise intensity is large, we will see again as in previous 
example that essential nonlinearity of the potential is suppressed by strong fluctuations and the 
evolution of the mean coordinate becomes exponential (see Fig. 4, plotted in the logarithmic scale, 
a = 1, b = 2, x = 2.5). 




(17) 
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V. CONCLUSION 



In the frame of this paper we have considered time evolution of averages of an overdamped Brow- 
nian motion. Extending the approach by Malakhov we have obtained exact characteristic scales 
of time evolution of averages (valid for arbitrary potentials and noise intensity). Some particular 
examples of time evolution of mean coordinate of overdamped Brownian motion in monostable and 
bistable potentials have been analyzed. 

It has been demonstrated that for purely noise-induced transitions (transitions over potential 
barriers) the time evolution of mean coordinate is a simple exponent with a good precision even for 
the case when the potential barrier height is comparable or smaller than the noise intensity. 

Also, there has been observed the effect of "linearization" of nonlinear system by a strong noise in 
the sense that time evolution of mean coordinate becomes purely exponential with increase of noise 
intensity. 
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FIG. 1. Evolution of the mean coordinate in the potential <P(x) — a; 4 /4 for different values of noise intensity. 
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FIG. 2. Evolution of the mean coordinate in the potential <P(x) — a; 4 /4 for different values of noise intensity. 
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FIG. 3. Evolution of the mean coordinate in the potential ${x) = x 4 /4 — x 2 for different values of noise intensity with initial 
distribution located in a potential minimum. 
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FIG. 4. Evolution of the mean coordinate in the potential ${x) = x 4 /4 — x 2 for different values of noise intensity with initial 
distribution located far from a potential minimum. 
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